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ABSTRACT 

We report on the detection of infrared light from the super-Earth 55 Cnc e, based on four occultations 
obtained with Warm Spitzer at 4.5 /im. Our data analysis consists of a two-part process. In a first 
step, we perform individual analyses of each dataset and compare several baseline models to optimally 
account for the systematics affecting each lightcurve. We apply independent photometric correction 
techniques, including polynomial detrending and pixel-mapping, that yield consistent results at the 
l-cr level. In a second step, we perform a global MCMC analysis including all four datasets, that 
yields an occultation depth of 131 ± 28 ppm, translating to a brightness temperature of 2360 ± 300 K 
in the IRAC-4.5/zm channel. This occultation depth suggests a low Bond albedo coupled to an 
inefficient heat transport from the planetary dayside to the nightside, or else possibly that the 4.5 /zm 
observations probe atmospheric layers that are hotter than the maximum equilibrium temperature 
(i.e., a thermal inversion layer or a deep hot layer). The measured occultation phase and duration 
are consistent with a circular orbit and improves the 3-cr upper limit on 55 Cnc e's orbital eccentricity 
from 0.25 to 0.06. 

Subject headings: planetary systems - stars: individual (55 Cnc, HD 75732) - techniques: photometric 



1. INTRODUCTION 

The nearby sixth magnitude naked-eye star 55 Cnc 
is among the richest exoplanet systems known 
so far, wit h five pl anetary co mpanions d e tected 



since 1996 (iButler et al. 19 93; iMarcv et al.[ 
McArthur et al. 2004; Wisdomll2005l: iFischer et all 



2002 



2008 



Dawson fc Fabrvckv .2010f ). The recent discovery of the 



of Mp = 
Rp = 2.17±0.10i?, 



transiting nature of 55 Cnc e, made independently in 
the visible with the MOST satellite (|Winn et all 120 lit ) 
and in the infrared with the Spitzer Space Telescope 
(jDemorv et al.lf201lL hereafter Dll), set this super-Earth 
among the most promising low-mass planets for follow- 
up characterization. 

A recent data reanalysis combining fifteen days of 
MOST monitoring and two Spitzer transit observations, 
allow ed us to refine 55 Cnc e's properties (jGillon et al.l 
120121 hereafter G12) that result in a planetary mass 
7.81 ± 0.56 M0 and a planetary radius of 
Although they do not uniquely con- 
strain the planetary composition, the mass and radius 
(and hence density) are consistent with a solid planet 
without a large envelope of hydrogen or hydrogen and 
helium. The planet could be a rocky core with a thin 
envelope of light gases. Another possible mass/radius 
interpretation could be an envelope of supercritical wa- 
ter above a rocky nucleus, where the exact amount of 
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volatiles would depend on the composition of the nucleus 
(G12). In this scenario, 55 Cnce would be a water world 
similar to the core of Uranus and Neptune. 

The super-Earth size of 55 Cnc e together with an ex- 
tremely high equilibrium temperature ranging between 
1940 and 2480 K, motivated us to apply for Spitzer Di- 
rector's Discretionary Time to search for the occultation 
of 55 Cnce at 4.5 /im. The goal was two-fold. First, the 
occultation depth provides an estimate of the brightness 
temperature, constraining both the Bond albedo and the 
heat transport e fficiency between the planetary day and 
night-side (e.g., iCowan fc Agoll[20Tll) . Second, a mea- 
surement of the occultation phase and duration provides 
a constraint on a po tential non-zero orbital eccentricity 
(jCharbonneaul l2003l ) that could be maintained by the 
interactions with the four other planets of the system. 

We present in this Letter the first detection of light 
from a super-Earth. The new Warm Spitzer 55 Cnc 
observations and corresponding data reduction are pre- 
sented in Section 2, while the photometric time-series 
analysis is detailed in Section 3. We discuss in Section 4 
the implications for our understanding of this planet, 
especially regarding the constraints on its atmospheric 
properties and orbital evolution. 



2. OBSERVATIONS 

Four occultation windows were monitore d by Spitzer 
in th e 4.5 /im channel of its IRAC camera ()Fazio et al.l 
[200l in January 2012. Table 1 presents the description 
of each Astronomical Observation Request (AOR). For 
each occultation, 6230 sets of 64 subarray images were 
acquired with an individual exposure time of 0.01s. All 
the data were calibrated by the Spitzer pipeline version 
S19.1.0 which produced the basic calibrated data (BCD) 
necessary to our reduction. For all runs except the first 
one, the new Pointing Calibration and Reference Sen- 
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sor (PCRS) peak-up mod^ was enabled. Because of the 
intrapixel sensitivity variability of the IRAC InSb detec- 
tors coupled to the point response functiorQ (PRF) un- 
dersampling, the measured flux of a point source shows 
a strong correlation with the intrapixel position of the 
star's center. This well documented "pixel-phase ef- 
fect" creates a correlated noise that is the main lim- 
itation on the phot ometric precision of Warm Spitzer 
(|Ballardet all l20Tol ). The new PCRS peak-up mode 
aims at mitigating this correlated noise by improving the 
telescope pointing's accuracy. 

We apply the same reduction procedure for all AORs. 
We first convert fluxes from the Spitzer units of specific 
intensity (MJy/sr) to photon counts. We then perform 
aperture photometry on each subarray image using the 
APER routine from the IDL Astronomy User's Librarj0. 
We compute the stellar fluxes in aperture radii ranging 
between 2.0 and 4.0 pixels, the best results being ob- 
tained with an aperture radius of 3 pixels for the first 
and fourth AOR, 2.8 pixels for the second AOR and 2.9 
pixels for the third AOR. We use background annuli ex- 
tending from 11 to 15.5 pixels from the PRF center. The 
center and full width at half maximum (FWHM, along x 
and y axes) of the PRF is measured by fitting a Gaussian 
profil e on each irnage u sing the MPCURVEFIT proce- 
dure ()Markwardtll2009[ ). 

Examination of the background time-series reveals an 
"explosion" of the flux in the third and fourth AORs, 
with a similar shape to the one observed in Dll. These 
two background "explosions" seem to be correlated with 
an increase of the measured FWHM along the x-axis. 
Because of the "pixel-phase" effect, changes in the PRF's 
shape create an additional correlated noise contribution 
(Gillon et al., in prep.). 

For each block of 64 subarray images, we discard the 
discrepant values for the measurements of flux, back- 
ground, x-y positions and FWHM using a lO-cr median 
clipping for the six parameters. We then average the re- 
sulting values, the photometric errors being taken as the 
uncertainties on the average flux measurements. At this 
stage, a 50-cr clipping moving average is used on the re- 
sulting light curve to discard totally discrepant subarray- 
averaged fluxes. The number of frames kept for each 
AOR are shown in Table 1. 

Figure 1 shows the raw lightcurve, the background flux, 
the x-y centroid positions and FWHM for each AOR. 
The improved stability of the telescope pointing achieved 
by the new PCRS peak-up mode can be easily noticed for 
the X time-series but seems less sharp along the pixel's y 
axis. 

3. DATA ANALYSIS 

3.1. Independent analysis of each AOR 

The aim of this step is to perform an exhaustive model 
comparison to determine the optimal baseline model for 
each AOR. The baseline model accounts for the time- 
and position- dependent systematic effects relevant to 
the IRAC-4.5/xm observations (see Dll, G12 and ref- 
erences therein). For this purpose, we individually an- 



° http://irsa.ipac . caltech.edu/data/SPITZER/docs/irac/pcrs. 
We followed here the IRAC instrument handbook terminology: 
'http : //irsa. ipac . caltech.edu/data/SPITZER/docs/iracj 
http: //idlastro . gsfc .nasa.gov/contents .html 



alyze each AOR by employing our adaptative Markov- 
C hain Monte Carlo ( MCMC) implementation described 
in iGillon et al.l (|2010[ ). We set the occultation depth as 
a jump parameter and fix the orbital period P, transit 
duration VF, time of minimum light Tq and impact pa- 
rameter b — a cos i/ Ri, to the values obtained from the 
global analysis of the system (G12). We further impose a 
Gaussian prior on the orbital eccentricity (e = 0.06±0.05, 
Dll), allowing the eclipse phase and duration to float in 
the MCMC. For each model, we run two chains of 10^ 
steps each. Throughout this work, we assess the conver- 
gence and good mi xing of the Marko v chain s using the 
statistical test from lGelman fc RubinI (|1992[ ). 

We first assume a baseline model based on a classical 
second order x-y position polynomial (Dll, eq. 1) to cor- 
rect the "pixel-phase" effect, added to a time-dependent 
linear trend. We then increase this basic baseline com- 
plexity by trying combinations of up to third- and fourth- 
order x-y position polynomials, second-order logarithmic 
ramp models and second-order time-dependent polyno- 
mials. The baseline models tested are therefore char- 
acterized by 7 to 20 free parameters, well constrained 
by the ~6200 measurements of each AOR. We finally 
compute the Bayesi an information criterion (BIG, e.g., 
IGelman et al.l I2003D for all combinations and choose, 
from the MCMC output, the baseline model that yields 
the highest marginal likelihood. We show the resulting 
individual lightcurves on Fig. 2 and the selected model 
along with the occultation depth for each AOR in Ta- 
ble 1. 

The correlated nois e affecting each ligh tcurve is taken 
into account following lGillon et al.l ()2010| ). A scaling fac- 
tor j3 from the comparison of the standard deviation of 
binned and unbinned flux residuals is determined during 
a preliminary MCMC run. This factor is then applied to 
the individual uncertainties of the flux time-series. The 
/3 values for each AOR are shown in Table 1. To obtain 
an additional estimation of the residual correlated noise, 
we conduct a residual-permutation bootstrap analysis on 
each lightcurve corrected from the baseline model, simi- 
lar to Dll. This part of the analysis yields parameters in 
good agreement with the results from our MCMC analy- 
ses, albeit with significantly smaller error bars, suggest- 
ing that the error budget is dominated by the uncertain- 
ties on the coefficients of the complex baseline model and 
not by the residual correlated noise contribution. 

The examination of the individual lightcurves obtained 
with more complex models lead to similar results to the 
lightcurves shown on Fig. 2, that are obtained with the 
baseline models selected in the previous step (Table 1). 
The resulting occultation depths are compatible within 
Icr, securing both our detection and the insensitivity of 
the occultation signal to our adopted method for system- 
atics correction. 

As in Dll and G12, we perfor m a Lomb-Scargle peri- 
odogram analysis ()Scarglelll982| ) on the residuals of our 
selected model for each AOR. Peaks with marginal sig- 
nificance are found at 69 min in the second AOR and at 
51min in the third AOR, while neither of the first nor 
fourth AOR's residuals reveal a periodic signal. We in- 
obscytSfiithese sinusoidal modulations in the baseline models 
of our second and third AOR and perform a new MCMC 
that yields a higher BIG value than our model selected 
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Fig. 1. — Flux with best-fit baseline models superimposed (top left), background (top right), centroid position on the x-axis (middle left) 
and y-axis (middle right), PRF FWHM on the x-axis (bottom left) and y-axis (bottom right), from first (top) to fourth AOR(bottom). 
The time-series are arbitrarily shifted vertically for clarity. 
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during the previous step of the analysis. We therefore 
neglect the sinusoidal term in both AORs. No hint of a 
~ 50min and ~ 100 ppm amplitude modulation similar 
to the one reported in Dll and G12 is found in the four 
datasets. 

We further build a "pixel map" to characterize the in- 
trapixel variability on a fine grid. This approach has 
been already de i nonstr ated bv lBallard et al.l ()2010D and 
iStevenson et al.l (|201lD as an efRcient method to remove 
the flux modulation due to the "pixel-phase" effect. The 
improved tracking accuracy brought by the new PCRS 
peak-up mode is indeed expected to increase the reso- 
lution of the pixel map, hence motivating this step of 
the analysis for our observations. More specifically, our 
implementation divides the area covered by the PRF in 
a grid made of 30x30 boxes and counts the number of 
out-of-eclipse datapoints that fall in a each box. If a 
given box has at least four datapoints spanning at least 
50% of the AOR duration, the individual fluxes are di- 
vided by their mean, otherwise the corresponding mea- 
surements are rejected. As this method could average out 
an eclipse signal located in the designated out-of-eclipse 
parts of the lightcurve, we repeat this procedure by grad- 
ually shifting the out-of-eclipse datapoints in time. On 
average, as compared to the first AOR (obtained with- 
out the PCRS peak- up mode), the measurements sample 
24% more time per box in the second AOR, 57% in the 
third AOR and 48% in the fourth AOR. The occultation 
depth and duration obtained using the pixel-map cor- 
rected lightcurves are in good agreement with the ones 
employing a polynomial baseline model. 

The different analysis techniques applied above yield 
detrended time-series in which the occultation signal is 
visible by eye with consistent phase and duration in three 
out of the four AORs (Fig. 2). The second AOR is the 
only one in which the detection is marginal (68±52ppm). 
The different approaches used to correct the photometry 
from the "pixel-phase" effect show that the intrapixel 
sensitivity does not explain this discrepancy. We notice 
that the (3 factor for this second AOR is ^-^50% larger 
than for the other AORs (see Table 1), suggesting a larger 
amount of correlated noise of instrumental or astrophysi- 
cal origin. Examination of the onboard temperature sen- 
sors readings and other external parameters in the FITS 
files do not reveal any unusual patterns. While an actual 
variability of the planet's emission cannot be ruled out, 
the marginal disagreement (1.2cr) of the second occulta- 
tion's amplitude relative to the three others is probably 
caused by this larger correlated noise. 

3.2. Global analysis 

The final step of our analysis consists of performing 
a global analysis, including all four lightcurves in the 
same MCMC framework, to constrain both the occulta- 
tion depth and orbital eccentricity. For this purpose, we 
assume Gaussian priors on 6, W, Tq, ^/ecosoj, yesino; 
and the stellar parameters based on t he posterior dis- 
tributions derived in Dll, G12 and in Ivon Braun et al.l 
(poll). We first run two MCMC (one with the occul- 
tation model and one without) to assess the robustness 
of our detection, using as input the four raw lightcurves. 
We employ the baseline models selected in the previous 
step and shown in Table 1. These two runs are composed 
of three chains of lO'' steps each. 



We find an occultation depth of 131 ± 28 ppm, and 
an eccentricity e < 0.06 (S-cr upper limit). The odds 
ratio computed using the BIC between the two models 
(with vs. without occultation) is ~ 10** in favor of the 
occultation model. The thermal emission from 55 Cue e 
is therefore firmly detected. 

To test the robustness of the eccentricity signal, we 
then perform an identical MCMC run but with the as- 
sumption of a circular orbit. The odds ratio between the 
circular and non-circular orbits is ~ 10'^ in favor of the 
circular case. The resulting occultation depth in this case 
is 122 ± 21 ppm, in excellent agreement with the value 
obtained for the non-circular case. The phase-folded oc- 
cultation lightcurve with the best-fit circular model is 
shown on Figure 3. 

4. RESULTS AND DISCUSSION 
4.1. Planetary Properties 

From the measured secondary eclipse depth we obtain 
a brightness temperature estimate of 2360 ± 300 K, us- 
ing a stellar blackbody emission spectrum, with T^ff — 
5196 ± 24K (von Braun ct al. 2011). At face value, the 
high brightness temperature suggests that cither 55 Cue e 
has both a low Bond albedo and an inefficient heat trans- 
portation from the day side to the night side, or the ob- 
servations in the IRAC 4.5 /zm bandpass probe layers in 
the atmosphere that are at temperatures higher than the 
equilibrium temperature (Figure 4) . Either scenario may 
explain why the brightness temperature is observed to be 
higher than the zero Bond albedo equilibrium tempera- 
ture, Teq — 1950 K, for a uniformly re-radiating planet. 

One interpretation could be that the planet is a rock 
with only a minimal atmosphere established through 
vaporization pchaefer fc Fcgicy 2009; Lcgcr et al. 200^ 
iCastan fc Menoull2blli ). Rocky objects in the Solar Sys- 
tem, e.g.. Mercury and the Moon, have low Bond albedos 
between 0.07 and 0.12. Lacking a thick atmosphere, a 
rocky super-Earth also does not have any efficient mean 
of transporting heat from the day-side to the night-side. 
We disfavor the bare rock scenario, however, because the 
55Cnce mass and radius measurements (Dll, G12) ex- 
clude a rocky composition similar to Mercury and Earth; 
to be a rocky planet with minimal atmosphere 55 Cue e 
would have to have the unlikely bulk composition of pure 
silicate. 

Alternatively, if 55 Cnc e has a substantial gas atmo- 
sphere or envelope, the scenario of inefficient heat redis- 
tribution is still supported by the T = 2360 K brightness 
temperature, as long as the radiation at 4.5 fim is not 
coming from a very hot layer, such as a deep layer or a 
thermal inversion layer. The interpretation of 55 Cnc e 
having a supercritical water envelope above a solid nu- 
cleus (Dll) fits with the inferred high temperature; the 
absence of clouds at such a high temperature results in a 
water envelope with a naturally low albedo. Although we 
cannot fully exclude the case of heat redistribution, for 
the water planet or even a gas atmosphere over a rocky 
core, the probe to deep hot layers or the thermal inver- 
sion would have to be extreme. A probe of extremely 
deep, hot atmosphere layers is unlikely because many 
of the gases likely to be present in the atmosphere, in 
particular CO2 and CO, have high absorption cross sec- 
tions in the spectral region of the IRAC bandpass be- 
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Fig. 2. — Left: lightcurves obtained during the independent analysis of each AOR (see Sect. 3) divided by the best-fit baseline models. 
Right: residuals for each AOR. AORs are displayed from first (top) to fourth AOR (bottom). Lightcurves are binned per 5 min. Individual 
detections in each AOR are consistent both in phase and duration. 
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Fig. 3. — Phase- folded occultation lightcurve resulting from our global analysis presented in Sect. 3. The best-fit circular model is 
superimposed in red. The lightcurve is binned per 8 min. 
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TABLE 1 

Individual AOR and global fit properties 



AORl 



AOR 2 



AOR 3 



AOR 4 



GLOBAL (prior on e) GLOBAL (e fixed) 



Observation date [UT] 
Observation window [UT] 
Number of measurements 
Baseline model* 
13 factor'' 

Occultation depth [ppm] 

\/ecosa) 

\/esinaj 



2012-01-18 
02:45 - 08:38 
6187 

1.25 
202 ± 54 

7+0.054 
-0.042 



0.037"' 



0.003 



+0.079 
-0.081 



2012-01-21 
01:16 - 07:09 
6190 

xy"^ + 1"^ 
1.89 

68± 52 

'-'•'^'^-0.108 

0.008+0;i|J, 



2012-01-23 
06:24 - 12:17 

6133 
xy"^ + t'^ +r'^ 
1.20 
107 ± 53 

..+0.086 
-0.077 
,+0.092 



0.036: 
-0.001"' 



0.100 



2012-01-31 
08:56 - 14:49 

6464 
xy^ +t'^ + 
4.27 
487 ± 56 
004+°-''*'i 
040+" '"^2 



434 lb 28 

U.U»'i_Q oj^g 

001+" "^* 



122 ±21 



* Baseline models are described by position (xy), time (t) and ramp (r) terms, with the order indicated in superscript. 
^ See Section 3 for details on the /3 scaling factor. 
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tween 4 and 5 /xm. A thermal inversion, while present 
in Solar System planet atmospheres and suggested to be 
prese nt in many similarly highly-irradiated hot Jupiters 
(e.g., iBurrows et 311120071 iKnutson et "aD l2008| ) . is also 
an extreme explanation because the temperature would 
have to be ~ 500 K above most of the rest of the atmo- 
sphere. In general, while we know the total amount of 
energy re-radiated by the planet, it is the interplay be- 
tween absorption of stellar radiation at short wavelength, 
the opacities at wavelengths at which the bulk radiation 
is emitted {^1-5 fim for a body at T~ 2000 K), and the 
opacity in our IRAC 4.5 fim. bandpass that sets the ob- 
served brightness temperature. 

Considering the observational uncertainty, the ob- 
served brightness temperature in the IRAC bandpass 
could be as low as 1830 K at the 2-a level. At this level 
of uncertainty, we can say that at least one or more of 
the following statements are true: 1) the planet's Bond 
albedo is low; 2) the planet has an inefficient heat trans- 
port from the day to night side; and/or 3) the IRAC 
4.5 ^m bandpass probes at atmospheric levels that are 
considerably hotter than the day-side equilibrium tem- 
perature. 

4.2. Orbital Eccentricity 

Our global MCMC analysis improves the 3-cr upper 
limit on 55Cnce's orbital eccentricity from 0.25 (Dll) 
to 0.06. However, even such a small eccentricity is un- 
likely because tidal interactions between the planet and 
star would probably damp an initial eccentricity that 
large in a few Myrs. Much larger initial eccentricities 
would be damped on similar timescales. However, high 
order dynamical interactions among the planets in the 
system may maintain an eccentricity for 55 Cnc e of a few 
parts per million, by comparison with theoretical stud- 
ies of multi-planet s ystems having close-in super-Earths 
(|Barnes et al.l |2010| ). Strong observational constraints 
on 55 Cnc e's eccentricity may even provide information 
about the planet's tidal response and internal structure 



(jBatygin et al.l [20091) . While a small eccentricity might 
have little direct influence on observation, the concomi- 
tant tidal dissipation within the planet can have dra- 
matic geophysical consequences, perhaps powering vig- 
orous volca nism and rcsupplying the planets atmosphere 
(| Jackson et al. 2008). The planet's proximity to its host 
star suggests, in fact, it may be shedding its atmosphere, 
and so active resupply may be necessary for long-term 
atmospheric retention. 

4.3. Future prospects 

From the first, landma rk detection of infrared 
light from a hot Jupiter (ICharbonneau et al.l 120051 : 
IDeming et al.|[20 05) and further occultation observations 
of more than t wo-dozen Neptune- and Jupiter-sized plan- 
ets (see, e.g., iDeming fc Seageil 12009) . Spitzer remains 
today the best instrument for exoplanet high-precision 
near-infrared photometry. At the dawn of the super- 
Earth sized planet discovery era (see, e.g.. IBatalha et ahl 
I20121 ). Spitzer 's continued legacy shows the need for 
keeping this observatory operational. With the future 
launch of the James Webb Space Telescope, the exo- 
planet community anticipates thermal emissions mea- 
sure ments for a numbe r of different super Earths (see, 
e.g.. lDeming et al.ll2009( ). hence improving our knowledge 
of this class of planets in a way similar to the achieve- 
ments made by Spitzer for hot Jupiters. 

We thank Julien de Wit, Renyu Hu, Sean Carey and 
Nikole Lewis for insightful discussions and the anony- 
mous referee for a report that improved the paper. This 
work is based in part on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy under a contract with NASA. We thank the Spitzer 
Science Center staff and especially Nancy Silbermann for 
the efficient scheduling of our observations. M. Gillon is 
FNRS Research Associate. 
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